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Abstract 

The design of hybrid peptide-solid interfaces for nanotechnological applications such as 
biomolecular nanoarrays requires a deep understanding of the basic mechanisms of peptide 
binding and assembly at solid substrates. Here we show by means of experimental and com- 
putational analyses that the adsorption properties of mutated synthetic peptides at semicon- 
ductors exhibit a clear sequence-dependent adhesion specificity Our simulations of a novel 
hybrid peptide-substrate model reveal the correspondence between proline mutation and bind- 
ing affinity to a clean silicon substrate. After synthesizing theoretically suggested amino-acid 
sequences with different binding behavior, we confirm the relevance of the selective mutations 
upon adhesion in our subsequent atomic force microscopy experiments. 

Introduction 

In the past few years, the interest in hybrid interfaces formed by "soft" molecular matter and "hard" 
solid substrates has rapidly grown as such systems promise to be relatively easily accessible can- 
didates for novel biosensors or electronic devices. The enormous progress in high-resolution mi- 
croscopy and in biochemical engineering of macromolecules is the major prerequisite for studies 
of hybrid systems and potential applications. ^'^ One particularly important problem is the self- 
assembly and adhesion of polymers, proteins, or protein-like synthetic peptides to solid materi- 
als such as, e.g., metals, ^'^ semiconductors,-"- carbon and carbon nanotubes,'^*^^ and silica.— ^^^ 



Peptide and substrate specific binding affinity is particularly relevant in pattern recognition pro- 
cesses.—^ Systematic experimental studies have been performed to investigate binding properties 
of individual amino acids in their binding behavior to selected materials. - Basic theoretical con- 
siderations of simplified polymer-substrate and protein-substrate models have predicted complex 
pseudophase diagrams.—^ 

In bacteriophage display experiments, only a few peptides out of a library of 10^ investigated 
sequences with 12 amino acid residues were found to possess a particularly strong propensity to 
adhere to (100) gallium-arsenide (GaAs) surfaces.^ The sequence- specificity of adsorption strength 
is a remarkable property, but the question remains how it is related to the individual molecular 
structure of the peptides. We expect that relevant mutations of sites in the amino-acid sequence 
can cause a change of the binding affinity. Indeed, one key aspect of our study is to show that 
proline is a potential candidate for switching the adsorption propensities to cleaned (100) silicon 
(Si) substrates. 

Silicon is one of the technologically most important semiconductors, as it serves, for example, 
as carrier substrate in microelectronics. For this reason, electronic and surface properties of Si are 
well investigated. This regards, for example, oxidation processes in air^^'^^ and water,— »^ as well 
as the formation of hydride surface structures and the Si-binding characteristics of small organic 
compounds. -^»2^ 

Results and Discussion 

To guide the design of peptide-silicon interfaces, we first performed extensive computer simula- 
tions of a novel hybrid model discussed below. For testing the theoretically revealed trends of 
adsorption propensity changes by selected mutation, we have synthesized the suggested specific 
mutants by means of multiple solid phase peptide synthesis. The theoretical predictions were 



subsequently verified in atomic-force microscope (AFM) experiments (see [[figure] [1][]1| and the 
detailed descriptions in the Supplementary Material). 



The hybrid model used in the computer simulations is composed of two parts contributing to 
the energy E(X) of a peptide conformation X: the energy of the peptide as represented by the 
implicit-solvent all-atom model introduced in Refs.— ^ and the interaction of the peptide with the 
substrate which is modeled in a coarse-grained way. The peptide model takes into account intrinsic 
excluded volume repulsions between all atoms, a local potential which represents the interaction 
among neighboring NH and CO partial charges, hydrogen bonding energy, and the interaction 
between hydrophobic side chains.— »2^ The substrate model consists only of atomic layers with 
surface specific atomic density and planar surface structure. In this simplified model, each pep- 
tide atom feels the mean field of the atomic substrate layers. The atomic density of these layers 
depends on the crystal orientation of the substrate at the surface. Based on these assumptions, a 
generic noncovalent Lennard- Jones approach for modeling the interaction between peptide atoms 
and surface layer is employed. ^'^^ We have studied this model by means of multicanonical com- 
puter simulations^^ which provide us with canonical statistics for any temperature T. The partition 
function is thus given by Z = / ^X^~^^^)/^^, where ^X is the formal functional integral measure 
for all possible conformations X in the space of the degrees of freedom. The statistical average 
of any quantity O is (O) = Z~^ / ^XO(X)^~^(^)/^^. In our simulations, the integral is estimated 
by an average over a large set of conformations (in each run about 10^ updates were performed) 
selected by multicanonical importance sampling. The precise modeling of the hybrid system and 
the multicanonical simulation methodology are described in the Supplementary Material. 

The peptide with the amino acid sequence SI [see [[figure] [2] []2[ a)] is a good example for 



the substrate specificity of adsorption. In recent comparative adsorption experiments, it could 
be shown that although SI binds strongly to GaAs(lOO), binding to Si(lOO) is very weak.^-^ In 
contrast, the adhesion is strongly increased, if the Si substrate is oxidized. This can clearly be 
[figure] [2] []2[ a), where AFM images of SI adhered at a de-oxidized (left) and an oxidized 



seen in 



Si(lOO) substrate (right) are shown. Peptide covered regions appear bright. A quantitative measure 
for the binding propensity is the peptide adhesion coefficient (PAC), which is the relative area of the 
surface covered by peptide clusters. ^'^ These PAC values are determined here by means of a cluster 



analysis of the respective AFM images. To reduce the dependence on the peptide concentration 
in solution, we here introduce the calibrated PAC (cPAC) as the ratio of PACs measured for the 
binding of the peptides to Si(lOO) and GaAs(lOO) substrates under identical conditions. GaAs(lOO) 
is chosen as a reference substrate, since the peptides considered here bind comparatively well to 
this substrate. The cPAC charts for SI in [[figure] [2] []2| ^b) clearly indicate the difference of binding 



affinity at cleaned and oxidized substrates. 

This is in contrast to sequence S3 [for sequence and AFM images see [[figure] [2] []2[ c)], which 



is a random permutation of SI with unchanged amino acid content. Surprisingly, the binding 
propensity of S3 to Si(lOO) was found to be much larger than that of SI.- In this case, the binding 
affinities at cleaned and oxidized Si(lOO) substrates are similarly strong, as the cPAC charts for S3 
[figure] [2] []2| ^b) show. In recent computational analyses of the solvent properties of these pep- 
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tides, we have shown that also the folding behaviors in solution exhibit noticeable differences. This 
is also true at room temperature, where in both cases the population of the structurally different 
native folds is rather small.— 

Another remarkable result of this former computational study is that the qualitative folding 
behaviors of SI and S3 are related to each other, if these sequences are pairwisely mutated at the 
position of proline, which occurs once in the sequences SI and S3.— The mutated sequence ST 
differs from SI only by the exchange of proline (P) at position 4 and threonine (T) at 9. Similarly, in 
S3', proline at 9 is exchanged with the aspartic acid (D) at 4, compared to S3 [see [[figure] [3] []3| ^a)]. 



These replacements were reasoned by the guess that the particular steric properties of proline, and 
thus its place in the sequence, influence the folding. It turned out that the folding behavior of ST 
in solution is close to that of S3, whereas S3' behaves rather like Sl.^^ Before we can address the 
question whether these results are also of essence for the adsorption behavior to Si(lOO), we need 
to discuss microscopic properties of de-oxidized Si(lOO) substrates. 

In our experiments, the Si(lOO) surfaces were first cleaned in a solution of ammonium fluoride 
(NH4F) and hydrofluoric acid (HF).-»^ Then, the peptide adsorption process took place in de- 
ionized water. This standard procedure ensures that the Si surface is virtually free of oxide and 



possesses strongly hydrophobic properties— »2^ (for sample preparation and experimental details 
see the Supplementary Material). The initial Si-F bonds after etching are replaced by Si-H bonds 
in the rinsing process in de-ionized water. After drying the sample, AFM scans of the surface 
were performed. Although the oxidation also proceeds in water,— »^ there are clear indications 
(maximum water droplet contact angle after removing the samples off the peptide solution) that 
the hydrophobicity of the Si samples remains largely intact during the peptide adsorption process. 
It is also known that Si surfaces are comparatively rough after HF treatment.— Thus, the reactivity 
of the surface is influenced by steps, which depend on the offcut and its directions. This renders 
an atomistic modeling intricate, even more as Si(100)-2xl surfaces are also known to form Si- 
Si dimers on top of the surface^ with highly reactive dangling bonds. From the considerations 
and the experimental preparations described above, it seems plausible that these bonds are mainly 
passivated by hydrogen, forming hydride layers. ^^»22i2^ it should be emphasized that under these 
conditions the surface structure of Si(lOO) is substantially different from oxidized Si(lOO) which is 
polar and in effect hydrophilic. ^^ An important result of [[figure] [2] []2| is that the binding of SI and 



S3 to oxidized GaAs(lOO) and Si(lOO) surfaces is virtually independent of the substrate type (cPAC 
'^ \). Thus, the top oxygen layer screens the substrate from the peptide. The different adhesion 
propensities to the clean (hydrated) substrates [see also [[figure] [2] []2| | lead to the conclusion that 



oxidation has not yet strongly progressed during the peptide adsorption process. We conclude that 
the key role of water is the slowing down of the oxidation process of the Si(lOO) surface, but for 
the actual binding process its influence is rather small. In particular, we do not expect that stable 
water layers form between adsorbate and substrate. 

These characteristic properties of HF treated Si(lOO) surfaces in de-ionized water effectively 
enter into the definition of our hybrid model of the peptide-silicon interface (for details see the 
"Model and Methods" section in the Supplementary Materials) which serves as the basis for our 
theoretical analysis and interpretation of the specificity of peptide adhesion on these interfaces. 

In order to quantify the degree of adsorption, we define the ratio of heavy (non-hydrogen) 
atoms located in a distance z/ < 5 A from the substrate, n^, and the total number of heavy atoms. 



Nk, as the adsorption parameter q = nk/Nk- The temperature dependence of its relative change 
by prohne mutation, A^(S^ -^ Sn') — {{q{Sn')) — {q{Sn))) / {q{Sn)) (with ^ = 1,3), is shown in 
[figure][3][]3[ b). The main result is that due to this selective mutation, the Si(lOO) adsorption 



affinity from SI to SI' increases [A^(S1 -^ SV) ^ 0.11 at T = 300 K], while it decreases by 
about the same amount as S3 is mutated to S3' [Aq{S3 -^ 83^) ^ -0.15 at T = 300 K]. This 
result is directly connected with the tendency to form secondary structures. In |[figure][3][]3[ c). 



the respective a-helix content [ratio of the dihedral Ramachandran angles of the inner 10 residues 
satisfying (j) G (-90°, -30'') and i/a g (-77°, -17'')] and j8-strand content [dihedral angles in the 
intervals ^ G (-150°, -90°) and Xf/ G (90°, 150°)] of the bound conformations are shown. We 
define a conformation as bound to the substrate, if at least 2% of the heavy atoms are within a 5 
A distance from the surface. There is a clear tendency that residues of SI and S3' are rather in a 
state and residues of S3 and SI' in j8 state. Yet the small secondary- structure contents are quite 
similar to what we found for the peptides in solution (without substrate), ^^ which were qualitatively 
consistent with analyses of CD spectra.- It is a noticeable result that here secondary structures are 
not stabilized near the cleaned Si(lOO) substrate, whereas in recent adsorption experiments of a 
synthetic peptide binding at silica nanoparticles, a stabilization of a-helices was observed. ^ - 
The experimental results shown in |[figure][3][]3| [d) confirm that the proline mutation of SI 



indeed increases the Si(lOO) binding affinity, while an analogous mutation decreases the binding 
strength of S3 by about the same value. The AFM images indicate the increased substrate coverage 
for SI ' and the decrease in the case of S3' [cp. with the corresponding images in [[figure] [2] []2| ^a) 



and (c), respectively]. By measuring the associated cPAC values, we find that AcPAC(Sl ^ Sl^) = 
cPAC(SlO -cPAC(Sl) ^ 0.27 and AcPAC(S3 -^ S30 ^ -0.25. This convincingly confirms our 
theoretical prediction from the hybrid-model simulations. 



Conclusions 

In summary, we have predicted by computer simulations and verified by AFM experiments that a 
selected proline mutation of short peptides facing a de-oxidized silicon substrate can substantially 
change the binding affinity in a very predictive and specific way. We could also show that this 
behavior is in part due to a qualitatively different folding behavior of the mutated sequences in 
the vicinity of the substrate. The proline position most likely also affects the aggregation prop- 
erties^ of the peptides and thereby indirectly again their binding characteristics. Building up on 
simulations of single-molecule behavior, as those discussed in the present manuscript, simulating 
coupled folding and aggregation while binding will therefore constitute a rewarding future project. 
Gaining deeper insights into the general principles of binding specificities is a first fundamental 
step towards the design of nanosensors with specific biomedical applications. Thus, the extension 
of our study to biomolecules is natural and the identification of unique bioprotein adsorption sig- 
nals in experiments with nanoarrays of several materials is a prerequisite for future applicability of 
such hybrid systems in biotechnology. 

The precise modeling of the hybrid system, the multicanonical simulation methodology and 
details of the peptide selection, the AFM experiments and the sample preparation are described in 
the Supplementary Material. 
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Figure 1 : Methodologies, (a) Sketch of the atomic-force microscope. The original AFM image 
exhibits SI peptide clusters on an oxidized lOx 10 [irn^ silicon substrate. The height of the highest 
cluster is 56 nm. (b) Computer simulations were performed with the BONSAI package that we de- 
veloped for Monte Carlo simulations of peptide models. The snapshot shows SI peptides forming 
helical segments near a silicon substrate. 
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Figure 2: Adsorption to cleaned and oxidized substrates. Exemplified AFM images of peptides 
(a) SI and (c) S3 adsorbed to cleaned (de-oxidized) and oxidized Si(lOO) surfaces. The AFM 
scale bar is 1 |Lim. (b) Calibrated peptide adhesion coefficients (cPAC) for SI and S3 adsorption 
to cleaned and oxidized Si(lOO) substrates. Single-letter code of amino acids occurring in the 
peptides of our study: A - alanine, D - aspartic acid, H - histidine, N - asparagine, P - proline, Q 
- glutamine, S - serine, T - threonine. 
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Figure 3: Reversed adsorption propensity of proline-mutated peptides, (a) Proline-mutated se- 
quences Sr and S3', (b) Adsorption parameter Aq as a function of temperature from our computer 
simulations. The prediction is that the prohne mutation of SI causes an increase of binding affinity, 
whereas mutating S3 leads to a decrease, i.e., proline mutation reverses the binding propensity of 
SI and S3, (c) a-helix content {na)h ^nd j8-strand content {np)\^ of bound conformations also ex- 
hibit a pairwise reversal of the tendency to form secondary structures. Exemplified conformations 
depicted in the insets are lowest-energy structures identified in the simulations (with rather small 
populations at room temperature) and represent the preferred trends in secondary- structure forma- 
tion: helical for SI and S3', sheet-like for SI' and S3, (d) Confirmation by AFM experiments at 
room temperature. The cPAC increases by proline-mutating SI and decreases by about the same 
value if S3 is mutated. AFM scale bar 1 |Lim. 
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Figure 4: Suggested TOC graphic. Designing novel hybrid peptide-solid interfaces for nanotech- 
nology requires a sound understanding of peptide binding and assembly at substrates. Here we 
show by means of experimental and computational analyses for a novel hybrid peptide-substrate 
system that the adsorption properties of mutated synthetic peptides at semiconductors exhibit a 
clear sequence-dependent adhesion specificity, which is mainly governed by the positions of cru- 
cial amino acids. 
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